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ABSTRACT 

We use the Cambridge stellar evolution code STARS to model the evolution 
of 5M and 7 Mq zero-metallicity stars. With enhanced resolution at the 
hydrogen and helium burning shell in the AGB phases, we are able to model 
the entire thermally pulsing asymptotic giant branch (TP-AGB) phase. The 
helium luminosities of the thermal pulses are significantly lower than in higher 
metallicity stars so there is no third dredge-up. The envelope is enriched in 
nitrogen by hot-bottom burning of carbon that was previously mixed in during 
second dredge-up. There is no s-process enrichment owing to the lack of third 
dredge up. The thermal pulses grow weaker as the core mass increases and 
they eventually cease. From then on the star enters a quiescent burning phase 
which lasts until carbon ignites at the centre of the star when the CO core 
mass is 1.36 Mq. With such a high degeneracy and a core mass so close to the 
Chandrasekhar mass, we expect these stars to explode as type 1.5 supernovae, 
very similar to Type la supernovae but inside a hydrogen rich envelope. 

Key words: stars: abundances, stars: AGB and post-AGB, stars: evolution, 
supernovae: general 



1 INTRODUCTION 

The primordial generation of stars, commonly referred to as Population-Ill stars or zero- 
metallicity stars, should have the composition of the interstellar (ISM) just after Big Bang 
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nucleosynthesis and hence have a negligible abundance of metals. It has been a popular 
belief that, in the absence of heavy elements and dust grains, cooling mechanisms are inef- 
ficient and favour the formation of massive or very massive stars. The s tandard mechanism, 
that accretion is terminated by radiation pressure on dust grains (e.g. 



Wolfire fc Cassinelli 



19871 ) in metal-rich gas, is not effective for gas with a primordial composition because there 



is no dust. Recently, it has been specul ated that accretion could instead be turned off 



through the formation of an H-II region (lOmukai fc Inutsuka 



tion pressure exerted by trapped 



( iPalla. Salpeter fc Stahler 



1983 



Lyq photons 



Yoshii fc Saio 



(Tan 



fc McKee 



2OO2I) or through the radia- 



20041 ) . It has also been shown 



19861 ) that even a small fraction of molecu- 



lar hydrogen can provide a significant contribution to cooling via rotational and vibrational 
transitions. The resulting Jeans mass of a pure H and He cloud could then be relatively small 
and may even fall below 0.1 M . Owing to this complexity and our lack of understanding 
of star formation, the initial mass function (IMF) of zero-met allici ty stars remains uncer 



Nakamura fc Umemura 



tain. U sing one- and two-dimensional hydrodynamical simulations, 
(]200ll ) showed that, depending on the initial density of the filamentary primordial gas cloud, 
there is an alternative result of the fragmentation of primordial filaments. For high density 
gas clouds, because the H 2 cooling is more effective owing to three-body reactions, filaments 
can contract and the fragmentation mass can be lo wered to 1M^. Hence, the y suggested a 
bimodal IMF with peaks close to 1 M & and 100 M . iJohnson &: Bromml (120061 ) have recently 
suggested that the formation of primordial low- and intermediate-mass stars is viable. It is 
likely that extremely low- to zero-met allicity AGB stars did form in the early Universe and 
so their evolution and contribution to the nucleosynthesis history should be investigated. 

In this paper, we describe the final stages of the evolution of 5M and 7 M & zero- 
metallicity stars, in particular the TP- AGB phase, and show that type 1.5 supernova are 
the probable fate of these stars. This shows that supernova explosions at zero-metallicity do 
not require initial stellar masses as high as in the case of solar metallicity and this can lead 
to important implications for galactic chemical evolution. 

The idea of a type 1.5 supernova or type l\ supernova is not new. Such supernova occur if 
a star's degenerate carbon/oxygen core grows up to near Ch andrasekhar m a ss before it loses 



Arnett 



(119691 ) . 



Ib en &: Renzini 



i ts envelope. This pos sibility was sugg ested, fo r insta nce, by 

(119831 ) . IWillsonl (120001 ) and recently by lZijlstral (120041 ) . As the core mass approaches 1.38 M & 
carbon ignites and the thermal runway in degenerate material cannot be delayed long enough 
to prevent an explosion from disrupting the entire star. Because the exploding star is a red 
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supergiant with a hydrogen-rich envelope, its spectrum and early light curve should closely 
resemble that of a supernova of type II. However, a substantial amount of radioactive Ni and 
Co is liberated by the exploding core, thus producing a late exp onential luminosity decline 



which could look like a type la supernova (jlben fc Renzinilll983l ). 



Lower metallicity stars have weaker stellar winds and thus their degenerate cores are able 
to grow up to near the Chandrasekhar mass and carbon ignition can lead to thermonuclear 
runaway and explosion. The mass-loss rate of low-metallicity AGB objects is uncertain but 
it is highly probable that it is lower than at sol ar metallicity. A faster core growth rate also 
increases the possibility of supernova Type 1.5. iGil-Pons. Gutierrez fc Garcfa-Berrol (120071 ) 
also reached the similar conclusion that supernovae of Type 1.5 are inevitable in the evolution 
of zero-metallicity stars between 5M Q and 7M Q , based on estimates of the mass-loss rate 
and core-growth rate. However, they did not compute the full evolution and estimated the 
two rates from the first few pulses. In this paper, we describe the full computation of 5 M & 
and 7 M Q models and show that carbon ignition at the degenerate core does occur. 



2 THE STARS CODE 



We use the Cambridge stellar evolution code stars to mode l the evolution of primordial 



intermediate-mas s stars. It was or i ginally written by 



Eggletonl (ll97ll ) and has been updated 



by many authors (IPols et alll995IUEldridge fc Toutl (20041 ) updated the opacity tables to use 
the latest OPAL calculations ( Iglesias fc Roger J 1996 ) and these new tables also account for 



changes in opacity with variations in the carbon and oxygen abundances. We use their zero- 
m etallicity opacities in this work (Eldridge, private communication based upon calculations 



by 



Ferguson et al 



20051 ). Unlike most codes which treat mixing in a separate step, this code 



solves the equations of stellar structure, nuclear burning and mixing simultaneously. 

Another unique feature of the STARS code is its use of a self-adaptive non-Lagrangian, 
non-Eulerian mesh. The mesh adapts so that mesh points concentrate in the physically 
important regions, such as burning shells and ionization zones, where things are changing 
most rapidly. During the asymptotic giant branch (AGB) phase, more mesh points are 
needed in the hydrogen and helium burning regions in order to resolve interaction between 
the two shells. Failure to resolve these regions properly can r e sult i n non-convergence of 



Straniero et al 



(119951 ) for a description of 



models or erroneous results. For example, see 
what happens with insufficient resolution. We use the AGB mesh spacing function described 
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by IStancliffe. Tout fc Pola (120041 ) to resolve thermal pulses. In the models presented here 
999 meshpoints were used. Such a large number of mesh points is needed to avoid resolution 
problems, not only during the AGB phase but also at the end of helium bu r ning and to avoid 

(boOO ). We do not 



Lau fc Tout 



the numerical problem for zero-metallicity stars described by 
include convective overshooting at any stage of the evolution. The inclusion of convective 
overshooting could lead to a larger core mass at the end of core helium burning but is unlikely 
to qualitatively chang e the evolution, though the surface abundances of metals could increase 



(IGil-Pons et al. 



2007) . We assu me there is no mass loss from the star. The mixing-length 



parameter (jBohm- Vitensdl 1 9581 ) a is 1.925 based on calibration to a solar model. The helium 
mass fraction is chosen to be 0.25 to reflect the prediction of primordial helium abundance 
from the observed deuteri um abundance, baryon den sity and a spectroscopic sample of 
extragalactic H-II regions (IFukugita fc Kawasakill2006l ). 



3 EARLY EVOLUTIONARY PHASES OF THE 7M MODEL 

The evolution of our 7 M & zero-metallicity star differs significantly from higher metallicity 
stars because of the absence of carbon, nitrogen and oxygen. Hydrogen cannot be burned 
through the CNO cycle, so it is burned via the proton-proton chain only. This chain is 
much less temperature dependent, so zero-metallicity stars are considerably hotter than 
their higher metallicity counterparts and their main-sequence lifetimes are much shorter. At 
the start of the main sequence, a convective core is driven by the pp-chain. The core ceases 
to be convective while hydrogen is still abundant. The temperature rises and becomes hot 
enough that carbon is produced by the triple-a reaction before hydrogen is exhausted. 
Because carbon is present in a hydrogen-rich region, the CNO cycle can now take place and 
dri ve a convectiye regi on in the core again. 



Chieffi et al. 



(120011 ) describe the central H and He burning of 4 — 7M stars. Our models 
agree with the characteristics they find. The convective core is much smaller than in metal- 
rich stars of similar mass. In our models it vanishes when the central mass fraction of 
H is 0.54, compared to 0.5 in their models. Our pp-chain-driven convective core is slightly 
smaller then theirs, while our CNO-cycle-driven convective core is slightly bigger than theirs. 
Secondly, as seen in Figure [TJ the He abundance increase s noticeably in 70 — 80 percent of 



the star by mass, in agreement with 



Chieffi et al. 



(120011 ). This is because the pp-chain is 



less temperature dependent than the CNO cycle. The helium abundance doubles out to 
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about 2M Q . The temperature gradient in the star is shallower than it would be in a higher 
metallicity star of the same mass. Consequently, hydrogen is burned in a more extended 
reg ion, so the helium abundance in creases over a larger part of the star. 

Siess. Livio fc Lattanziol (120021 ) describe the AGB phase of a 7 M & zero-metallicity star 



too. We compare our core masses a nd su rface abundances at this early AGB phase with 



theirs and the work of 



ChiefE et al. 



(120011 ) in Table 1. Our models are richer in nitrogen by 
a factor of ten and less abundant in carbon and oxygen by about a factor of ten and thirty 
respectively than those of Siess et al (2002). This suggests that their second dredge up is 
much deeper than ours and may even dip into the helium shell. This would also explain the 
drop in surface helium abundances for their 7M fl model c ompared to their own 5 M & model. 



Our results are closer to those of 



ChiefE et al 



(1200 ll ). The carbon and oxygen surface 
abundances agree very well for the 7 M Q models. One signiScant difference is that our 7 M Q 
model has a surface abundance of nitrogen about ten times higher. Comparison with their 
mass fraction proSles reveals that nitrogen is produced in a more extended mass range in 
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Figure 2. Mass fraction profile of the 7 Mq model at the end of central He burning. Note the significant amount of carbon 
and nitrogen produced at the hydrogen burning shell. 



Physical quantity This work Siess et al. (2002) Chieffi et al. (2001) 



core mass/AfQ 


1.0301 


1.0268 


0.9875 


4 He 


0.3821 


0.3764 


0.369 


12 C 


2.64 x 10~ 6 


2.41 x 10~ s 


2.08 x 10~ 6 


14 N 


2.28 x 10" 8 


1.44 x 10~ 9 


1.59 x 10~ 9 


16 Q 


3.85 x 10" 9 


9.39 x 10~ 8 


2.88 x 10~ 9 



Table 1. Comparison of core masses and surface abundances by mass fractions of the 7 Mq zero-metallicity early AGB models. 

our model, (compare Fig. [2] with Fig. 7 in their paper) so more of it is made and hence 
dredged up to the surface. This is because the hydrogen shell is thicker and the CNO cycle 
proceeds over a wider range of mass in our models. 

4 THE LATE AGB PHASE OF THE 7 M MODEL 

We have continued the evolution of the 7 M model through its entire thermally pulsing (TP) 
phase without any mass loss. There are 590 thermal pulses in 1.1 x 10 5 yr. The interpulse 
period is about 700 yr for the first few pulses and decreases to about 100 yr at the end of 
the TP phase. The core mass is 1.04 M & when the first thermal pulse starts. Our interpulse 
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Gil-Pons et al. 



(120071 ). Also, our 



periods are significantly shorter than the 2,700 yr found by 
core mass is bigger than their model without overshooting by 0.08 M@. In fact their model 
with overshooting has an interpulse period and core mass much closer to ours. The maximum 



helium luminosity of our model never exceeds 10 5 5 Lm. The pulses are to o weak to lead to any 



third dredge-up, in agreement with the work of 



Gil-Pons et al. 



( 120071 ). Carbon previously 



brought to the s urface dur ing second dredge-up has been converted to nitrogen by hot- 



bottom burning (llbenlll975l ). At this point, the surface carbon and nitrogen abundances are 
5.5 x 10~ 7 and 2.7 x 10^ 6 by mass. In the absence of third dredge-up, hot-bottom burning has 
reduced the surface carbon abundance during the TP-AGB nor have any s-process elements 
been brought up to the surface. 

The issue of third dredge-up in AGB stars has been a contentious one for some time. 
It has been postulated that thi rd dredge-up only happens in stars above a certain critical 



met a 



( 12002 



licity (iKomiya et al.l 120071 ) . While our models ag ree with the work o f 



Gil- 



(2007), who find that third dredge- up is absent, both 



ChiefS et al. 



(120011 ) and 



Pons et al 



Sicss et al 



) find that th i rd dre dge up does occur in their zero-metallicity intermediate mass AGB 



stars. 



Chieffi et al 



Her wig et al 



200 ll ) treated the convective boundaries according to the prescription of 



19971 ). who use a mixing scheme so efficient that the composition discontinu- 
ity between the two burning shells is smoothed out. This seems to indicate the efficiency of 
third dredge-up depends on the treatment of convectio n and the inclu s ion o f extra-mixing 



Gil-Pons et al 



(120071 ) find that the 



mechanisms such as convective overshooting. However, 
total amount of mass dredged up is very small even when overshooting is included. Prescrip- 
tions that overshoot into the processed core generally raise the metallicity and make the 
behaviour more like that of stars of higher metallicity that do undergo deep third dredge up 



( iStancliffe. Tout fc Pols 



els of super- AGB stars, 



2004: 



Stancliffe. Izzard fc Toutll2005l ). For example, in recent mod- 
D( . 1 1 ertv fc Lattanziol ( 200fil ) find that a 9.5 M Q star has a dredge- 



up efficiency A of 0.7. Dredge up efficiency is defined by the amount of H-exhausted core 
matter mixed int o the envelope divide d by the amount of core growth during interpulse 
period. However, ISiess fc Pumol (120061 ) find no third dredge up. As they highlighted, the 
occurrence (or not) of third dredge-up depends sensitively on how one treats the convective 
boundaries, as well as whether one includes additional mixing mechanisms. The stars code 



uses an arithmetic scheme for determining the diffusion coe fficient for the mixing. 



typically given deeper dredge- up than found in other codes (iStancliffe. Tout fc Pols 



t has 



2004 : 



Stancliffe. Izzard fc Toutll2005l ). so it is significant that we find no dredge- up in these models, 
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when others do. It may be that we do not find TDUP because we do not apply extensive 
extra-mixing and we shall investigate this in future work. 

Because helium burning already proceeds at a relatively high rate in the hotter burning 
shells during the interpulse period, the jump in helium luminosity during pulses is small 
compared to their higher metallicity counterparts. Later pulses are weaker and they cease 
altogether when the core mass reaches 1.1 M (see Figure [3]). The star then enters a quiescent 
evolutionary phase for about 1.8 x 10 5 yr while the hydrogen and helium burning shells grow 
outward without any thermal pulses. The sta r reach es carbon ignition after 2.9 x 10 5 yr, 



much faster than estimated by 



Gil- Pons et al. 



(2007^) at 1.2 x 10 6 yr. This is because their 



estimate was based on the core growth rate of the first few pulses. However, although our 
core growth rate is 4.7 x 10 -7 M yr -1 , which agrees with their rate for the first few pulses, 
it increases with time and is 2.0 x lO _6 M yr _1 when it explodes. 

We have compared a 7 M solar metallicity model with a 5 M zero-metallicity model 
with the same core mass in order to explain the weak thermal pulses of zero metallicity 
stars. We did not compare directly with the 7M zero-metallicity because its core mass is 
much larger than that of a 7 M solar metallicity model. We find that the zero-metallicity 
star has a much thinner intershell, helium-rich layer both in terms of mass and radius. 
There is less helium to be burnt during the thermal pulses, so they are much weaker and 
the interpulse period is much shorter. Because of the shorter interpulse period, the helium 
shell does not cool down as much as the higher metallicity model, so the temperature in the 
helium shell increases after each pulse. This is in contrast to the solar metallicity stars where 
the helium shell temperature during the interpulse goes down after each pulse. Eventually, 
the temperature becomes hot enough that helium burning can proceed smoothly without 
any pulses. This is because helium burning is much less temperature sensitive at higher 
temperatures. At T = lO 8 ^, the 3a reaction rate is pro portional to T 40 while at T = 2 x 10 8 K 
it is proportional to T 18 5 (cf. the stability criterion of lYoon. Langer van der Sluyal2004l ) . 



The hydrogen burning shell is much hotter and hence closer to the core in the zero-metallicity 
star. The higher temperature can cause earlier ignition of the pulse and hence the thinner 
intershell and helium burning shell. The hotter temperature can be attributed to the lower 
metallicity content of these stars, because of lower CNO abundances in the burning shell 
and lower opacity at the surface of the stars. 

Carbon ignition occurs at the centre under degenerate conditions when the core mass 
reaches 1.36 M & . The carbon luminosity rises rapidly and is soon followed by the breakdown 
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Figure 3. The end of the thermal pulses of the star for the 7Mq model. The pulses grow weaker and eventually stop as 
indicated by the variation of helium luminosity. Top: Model without mass loss. Bottom: Model with Reimers' mass loss. The 
difference between the two models is very small. 



of the evolution code because of the thermonuc lear runaway. We plot a carbon ignition curve 



as described by 



Martin. Tout fc Lesaffrd (120061 ) in Fig. [H Above and to the right of the solid 



line of the figure we have the right conditions for carbon ignition to drive a thermonuclear 
runaway. On the same axes, we plot the evolution of the internal temperature against density 
for the models leading up to ignition. Carbon ignites at the centre of the degenerate core 
because of the high density. Fig. H] shows that the core ignites carbon degenerately at the 
centre before any other part of the star. In particular, the burning shell is not hot enough. 
The subsequent thermonuclear runaway (similar to a type la supernova) releases sufficient 
energy to blow the whole star apart. Despite the envelope, we would expect the explosion 
mechanism of this star to be very similar to a Type la supernova, so we can estimate the 
nucleosynthetic yield from the explod i ng cor e. According to the deflagration model with zero 
metallicity (W70) of llwamoto et al.l (119991 ). the nucleosynthesis products of a supernovae 
type la would be 5.08 x 10" 2 M of 12 C, 3.31 x 1O~ 8 M of 14 N and 0.133 M of 16 0. As 
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described by lLau. Stancliffe fc Tout 



(120071 ). the yield of nitrogen from the envelope is of 
the order of 10~ 5 M® and for carbon 10 -6 M Q with an even lower oxygen yield, so we can 
conclude that the carbon and oxygen yields from the envelope are insignificant compared to 
the supernova yields. In the deflagration model nitrogen is mainly released by the envelope. 



However, if the explosion mechanism were delayed detonation the nitrogen yield o: 



Iwamoto et al. 



the core 



19991). In this 



would be about 2 x 1O~ 4 M (WDD1, WDD2 or WDD3 of 
case, the explosive yield would be higher than the yield from the connective envelop e and 
may be the source of nitrogen for the N-enhanced stars described by lSpite et al.l (120051 ) . The 
carbon yield drops to about 1O _2 M and the oxygen yield to about 7 x 1O _2 M . In both 
scenarios, the two biggest yields are silicon and iron. The 28 Si yield is 0.142 M to 0.272 M & 
and 56 Fe is 5.87 M & to 0.695 M & . Such stars are important for the iron contribution in the 
early Universe. The composition of material ejected by t he supernova can be very different to 
that released during binary interaction, as described by lLau et al.l (120071 1 when the envelope 
is lost before the ignition of carbon. 

In the above estimates, we have ignored the nucleosynthesis that may take place in 
the envelope during the explosion. Unlike a Type la supernovae which does not have a 
hydrogen-rich envelope, this star could produce extra nucleosynthesis during its explosion 
just as a Type II supernova does when a shock wave sweeps through the envelope. The 
presence of an envelope also makes fallback a possibility. The ejecta can then be enriched 
by a-capture isotopes and neutron-processed isotopes. If the star is in a binary system, r- 
processed elements formed during explosion can pollute the companion star. This ma y be a 
source of the double r/s-processed enriched halo stars, as suggested by IZijlstral (120041 ) . even 
though s-process elements are not brought to the surface in our models . 



5 MASS-LOSS RATE AND THE FATE OF THE 7M MODEL 



We have also evolved a 7M model with Reimers' mass loss. The evolution of this model 
is almost identical to that of the one without mass loss. Minor differences are that the 
thermal pulses stop about 1,000 yr earlier and the helium luminosity is generally lower, but 
by le ss than 0.1 dex, for the mo del with mass loss (see figure [3]). The Reimers' mass-loss 
rate (IKudritzki fc Reimers! 1 19 781 ) is 



M R = -4.0 x 10- 13 V iL/ ^!l /RQ) M Q yr-\ 
' (M/M ) OJ ' 



(1) 
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Figure 4. Interior of the 7.0 Mq model. The dotted line indicates where carbon burning produces energy faster than it is lost 
via neutrinos. We add the internal structure profile with the core at the highest density and the burning shell at the highest 
temperature. The dashed line is at a slightly earlier time and the solid line is at the end of the evolution. Carbon ignites 
degenerately at the centre before the shell is hot enough to ignite carbon. The long-dash line at the top left hand corner is the 
boundary between non-degenerate and degenerate regions. 



where L is the luminosity of the star, M its mass and R its radius. We uses 77 = 1, which is 
likely to be an overestimate of the mass loss in low-metallicity stars. The maximum mass- 
loss rate, at the end of evolution, when the luminosity and radius are highest, is 8.1 x 
lO _6 M yr _1 . The total mass lost from the star since the first thermal pulse is then about 
1.0 M & and the average mass-loss rate is 3.4 x 10~ 6 M ^yr~ 1 . We have also estim ated the 
mass- loss rate caused by the onset of radial pulsations (IVassiliadis fc Woodlll993l ). Because 
the star is much more compact than higher metallicity stars of similar mass, R/R & never 
exceeds 10 2 7 , so the maximum pulsation period is about 350 d. This p eriod is significantly 



Vassiliadis fc Wood 



below the period for the onset of the superwind phase discussed by 
(119931 ) . Therefore the mass- loss rate driven by pulsations is much lower than that of Reimers 
prescription. 
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The Reimers' rate was originally calibrated with M supergiants and so has no obser- 
vational support for AGB s tars s o we briefly conside r other mass-loss prescriptions such 



those of 



Schroder fc Cuntzl (120051) and 



The formula given by 

\ 3.5 



Schroder k Cuntzl d2005h is 



Blockerl Il995l ). which give higher mass- loss rates. 



M 



elf 



SC 



4,000K 



1 



9e 



4,300# 



Mi 



R- 



(2) 



where T e g is the effective temperature of the star, g is its surface gravity and g Q is the surface 
gravity of the Sun. This prescription gives a mass-loss rate that is 10.8 times that of the 
Reimers' prescription but still applies to non-pulsating giants and so may not be applicable 
to AGB stars. 



Blockerl (119951 ) gives the mass-loss rate as: 



M B = -4.83 x 10 



r2.7 



R • 



(3) 



For solar metallicity sta rs, this gives consistently larger rates during the TP-AGB phases 



.e.g. 



Gallart et al. 



the Magellanic Clouds, 



2005) so, based on the n umbers of luminous lithium rich AGB stars in 



Ventura et al. 



(120001 ) suggested the use of r] = 0.02 for the Reimers' 
mass-loss rate when using this prescription. With this modification, the mass-loss rate is 
about 15 times the standard Reimers' rate and the timescale for the loss of the envelope 
is only 1.1 x 10 5 years. Both these formulae give a mass-loss timescale shorter than the 
evolution time and the envelope could be lost slightly before carbon ignition. However, all 
our above estimates have neglected the eff e ct of metallicity. We can apply the commonly 



used scaling, suggested by 



Nugis &: Lamersl (120001 ) for hot stars, that 



0.5 



M(Z) = M(Z Q ) ( 



(4) 



where M(Z Q ) is the mass-loss rate for solar metallicity and Z is the surface metallicity. The 
scaling arises from the assumption that stellar winds are line driven. With lower surface opac- 
ity at lower metallicity there are weaker winds. However, while there is general agreement 
that mass lo ss falls with metal licity, there is a range of suggested values for the exponent. 
For example, 



Vink et al. 



3) 



(l2000f ) suggested M(Z) = M(Z )(^) 64 for B supergiants. Based 



o ■ 



on equation (4), as shown in table 2, all the mass-loss timescales ar e significantly long e r than 



the actual evolution time. Our estimated mass-loss timescale for 
is roughly 10 times shorter than the timescale given in 



Schroder fc Cuntzl (120051 ) 



Gil- Pons et al. 



(120071 ) . Our surface 



abundances are very close to their model, so the differences in timescale are because the star 
grows larger as it evolves. 
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The above scaling only applies to a radiation-driven wind and the relation comes from 
the fact that the efficiency of dust formation is reduced at low metallicity. However, as shown 
previously, the pulsation driven mass-loss rate is very low, even if the scaling does not apply 
Also, different mass-loss prescriptions may not be scalable by this relationship. Nevertheless, 
we plot the evolution time and mass-loss timescale with and without scaling in Figj6j The 
timescales for the star to lose its envelope are all significantly longer than the actual time 
the star takes to evolve up to the point of carbon ignition with the scaling. Even without 
scaling, the timescales are similar to, or slightly shorter than the evolution time, so even a 
weak scaling with a much higher surface metallicity means that the star does not lose its 
envelope before exploding. The limiting mass-loss rate is 1.9 x 10 5 M Q yr~ 1 . The unsealed 
Reimers' prescription gives a mass-loss rate much lower than the limiting rate. Unless the 
current prescriptions greatly underestimate the mass-loss rate from AGB stars, the fate of 
a 7M zero-metallicity star is to explode as type 1.5 supernova. However, whether a star at 
low metallicity can reach such a mass loss rate is not clear yet, due to lack of observational 
data. 



6 MASS-LOSS RATE AND THE FATE OF THE 5M© MODEL 

We have also evolved a 5M model with Reimers' mass loss with rj = 1. The evolution is 
very similar to that of the 7M star. When thermal pulses begin, the core mass is only 
0.92 M Q and pulses cease when the core has grown to 1.05 M . Like the 7M star, it also 
enters a quiescent phase until carbon ignites degenerately in the centre of the core when 
the core mass reaches 1.36 M (see Fig. [5]). The total time from the onset of thermal pulses 
to the onset of carbon ignition is 1.2 x 10 6 yr, much longer than the 7M star because the 
whole thermally-pulsing AGB phases lasts much longer. The surface metallicity of the 5 M & 
model is much lower because the convective envelope does not reach as deep during second 
dredge-up. The surface CNO abundances by mass fraction near explosion are 2.3 x 10~ 10 , 
5.1 x KT 9 and 3.2 x 10 -11 . 

The total mass lost from the star is 1.4 M , with an average mass-loss rate of 1.1 x 
10 6 M yr _1 . We have made a similar calculation of the mass-loss timescale based on different 
mass-loss rates (see Table 2). The limiting mass-loss rate is 3.0 x 10 6 M yr _1 so the possibility 
that the envelope is lost before carbon ignition is slightly higher because the core needs more 
time to grow. However, if we assume the mass-loss rate scales with metallicity, the mass-loss 
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Figure 5. Interior of the 5.0 Mq model. The dotted line indicates where carbon burning produces energy faster than it is 
lost via neutrinos. We add the internal structure profiles with the core at the highest density and the burning shell at highest 
temperature. The dashed line is at a slightly earlier time and the solid blue line is at the end of the evolution. Carbon ignites 
degenerately at the centre before the shell is hot enough to ignite carbon. The long-dash line at the top left hand corner is the 
boundary between non-degenerate and degenerate regions. 
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evolved time/yr 
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Time to lose envelope 
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5 
7 


No 
No 


1.2 X 10 6 
3.0 X 10 5 


3.1 X 10 6 
1.7 x 10 6 


3.7 X 10 5 
1.5 x 10 5 




4.9 x 10 5 
1.1 x 10 5 


5 
7 


Yes 
Yes 


1.2 x 10 6 
3.0 x 10 5 


6.0 x 10 9 
1.4 x 10 8 


6.9 x 10 8 
1.3 x 10 7 




9.0 x 10 s 
9.5 x 10 6 



Table 2. Timescales associated with loss of the envelope for different mass-loss rates compared with the actual evolution time. 
The top section has no metallicity-scaling of the mass loss is used. The bottom section uses the relation M(Z) = M(Z e )(-f-) - 5 . 



timescale is again much longer than the evolution time. As in the case of the 7M Q star, we 
can be fairly confident that carbon ignition and the following supernova do occur, unless 
the effect of overshooting or other extra mixing mechanisms such as rotation increase the 
surface metallicity and hence the mass-loss rate. 
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Figure 6. Timescales to lose the envelope using different mass-loss rates compared with the actual evolution time. The solid 
bars are when no metallicity scaling is used while the broken bars give the longer timescales involved when the mass loss scaling 
M(Z) = M(Zq)(-^—) 0,5 is applied, ft shows that only a small scaling is needed for the mass-loss timescales to be greater than 
the evolution time. 



7 CONCLUSION 

We have shown that the fate of high-mass AGB primordial stars is to ignite carbon de- 
generately at their centres and explode as supernovae with behaviour similar to a Type la 
but with a hydrogen-rich envelope because the mass-loss rate is low for these stars. Such 
supernovae enrich the early Universe with metals such as iron, nickel and carbon. Whether 
this occurs depends on the mass-loss rate from the star. A high mass-loss rate can cause the 
star to lose its envelope before carbon ignition and end its evolution as a white dwarf. For 
the 7 M Q star, the critical average mass-loss rate is 1.9 x 10 5 M Q yr _1 while for the 5M , 
it is 3.0 x 10 6 M Q yr _1 . So far there is no observational support for any of the proposed 
mass-loss rates at low metallicity. In order to be certain that these stars explode, we require 
the mass-loss rate of these stars to be less than about one-third of the solar rate. Because of 
the low surface metallicity of these objects, the mass-loss rate should be low enough unless 
the surface metallicity of the stars is substantially increased by extra mixing or the current 
mass-loss metallicity scaling is very wrong. 
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We have also shown that for extremely metal poor stars, the strength of thermal pulses 
is weak and there is a lack of third dredge-up. Eventually the thermal pulses disappear and 
the core growth rate is much faster. This is important for a supernova type 1.5 to occur 
because the core can grow much faster than the time it takes for star to lose its envelope. 
The lack of third dredge-up also has important implications for the contribution of AGB 
stars to the chemistry of the early Universe, particularly for s-process isotopes which we are 
not produced without it. Further models of the TP- AGB phases of very low metallicity stars 
are necessary to determine the lowest metallicity at which third dredge up begins. 

8 ACKNOWLEDGMENTS 

We thank Albert Zijlstra for his useful comments. HBL thanks PPARC for his Dorothy 
Hodgkin Scholarship. RJS and CAT thank Churchill College for their Fellowships. 

REFERENCES 

Arnett D. W., 1969, Ap&SS, 5, 180 
Bldcker T., 1995, A&A, 297, 727 

Bohm-Vitense E., 1958, Zeitschrift fur Astrophysik, 46, 108 

Chieffi A., Dommguez I., Limongi M., Straniero O., 2001, ApJ, 554, 1159 

Doherty C. L., Lattanzio J. C, 2006, Memorie della Societa Astronomica Italiana, 77, 828 

Eggleton P. P., 1971, MNRAS, 151, 351 

Eldridge J. J., Tout C. A., 2004, MNRAS, 348, 201 

Ferguson J. W. Alexander D. R., Allard F., Barman T., Bodnarik J. G., Hauschildt P. H., 

Heffner-Wong A., Tamanai A., 2005, ApJ, 623, 585 
Fukugita M., Kawasaki M., 2006, ApJ, 646, 691 
Gallart C, Zoccali M., Aparicio A., 2005, ARA&A, 43, 387 
Gil-Pons P., Gutierrez J., Garcia-Berro E., 2007, A&A, 464, 667 
Herwig F., Bloecker T., Schoenberner D., El Eid M., 1997, A&A, 324, L81 
Iben Jr. I., 1975, ApJ, 196, 525 
Iben Jr. I., Renzini A., 1983, ARA&A, 21, 271 
Iglesias C. A., Rogers F. J., 1996, ApJ, 464, 943 

Iwamoto K., Brachwitz F., Nomoto K., Kishimoto N., Umeda H., Hix W. R., Thielemann 
F.-K., 1999, ApJ, 125, 439 

© 0000 RAS, MNRAS 000.ITUT81 



The explosive end of Pop-Ill AGB stars 17 
Johnson J. L., Bromm V., 2006, MNRAS, 366, 247 

Komiya Y., Suda T., Minaguchi H., Shigeyama T., Aoki W., Fujimoto M. Y., 2007, ApJ, 
in press 

Kudritzki R. P., Reimers D., 1978, A&A, 70, 227 

Lau H. B., Stancliffe R. J., Tout C. A., 2007, MNRAS, 378, 563 

Lau H. B., Tout C. A., 2006, in ASP Conf. Ser. 353: Stellar Evolution at Low Metallicity: 

Mass Loss, Explosions, Cosmology p. 307 
Martin R. G., Tout C. A., Lesaffre P., 2006, MNRAS, 373, 263 
Nakamura F., Umemura M., 2001, ApJ, 548, 19 
Nugis T., Lamers H. J. G. L. M., 2000, A&A, 360, 227 
Omukai K., Inutsuka S.-i., 2002, MNRAS, 332, 59 
Palla F., Salpeter E. E., Stahler S. W., 1983, ApJ, 271, 632 
Pols O. R., Tout C. A., Eggleton P. P., Han Z., 1995, MNRAS, 274, 964 
Reimers D., 1975, Mem. Soc. R.Sci. Liege, Ser. 6, 8, 369 
Schroder K.-P, Cuntz M., 2005, ApJ, 630, L73 
Siess L., Livio M., Lattanzio J., 2002, ApJ, 570, 329 

Siess L., Pumo M. L., 2006, Memorie della Societa Astronomica Italiana, 77, 822 

Spite M., Cayrel R., Plez B., Hill V., Spite F., Depagne E., Frangois P., Bonifacio P., 

Barbuy B., Beers T., Andersen J., Molaro P., Nordstrom B., Primas F., 2005, A&A, 430, 

655 

Stancliffe R. J., Izzard R. G., Tout C. A., 2005, MNRAS, 356, LI 
Stancliffe R. J., Tout C. A., Pols O. R., 2004, MNRAS, 352, 984 

Straniero O., Gallino R., Busso M., Chieffi A., Raiteri C. M., Limongi M., Salaris M., 1995, 
ApJ, 440, L85 

Tan J. C, McKee C. F., 2004, in ASP Conf. Ser. 322: The Formation and Evolution of 

Massive Young Star Clusters p. 263 
Vassiliadis E., Wood P. R., 1993, ApJ, 413, 641 
Ventura P., D'Antona F., Mazzitelli I., 2000, A&A, 363, 605 
Vink J. S., de Koter A., Lamers H. J. G. L. M., 2000, A&A, 362, 295 
Willson L. A., 2000, ARA&A, 38, 573 
Wolfire M. G., Cassinelli J. P., 1987, ApJ, 319, 850 
Yoon S.-C, Langer N., van der Sluys M., 2004, A&A, 425, 207 
Yoshii Y., Saio H., 1986, ApJ, 301, 587 

© 0000 RAS, MNRAS OOO.ITUT81 



18 H.B. Lau, R.J. StancUffe & C.A. Tout 
Zijlstra A. A., 2004, MNRAS, 348, L23 



